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1 Introduction 

Quantum-Chromo-Dynamics (QCD) has become increasingly important to astro¬ 
physics during the last few decades. QCD is becoming especially ubiquitous in the 
study of compact objects, most specifically in the study of neutron stars (NSs). The 
hypothesis that quark matter (specifically, matter made of up, down, and strange 
quarks) could be the most stable state of matter in the universe [2j |3] has revealed 
the possibility for the existence of compact objects made entirely of deconhned quark 
matter namely, quark stars. The extremely dense cores of NSs could allow for the 
possibility of QCD phase transitions [I] leading to the formation of these quark stars. 

If the phase transition is an exothermic process, the energy could manifest itself 
observationally as an energetic event that could have drastic consequences to high- 
energy astrophysics. The exothermic release of energy by a QCD phase transition in a 
NS, resulting in a quark star, may lead to a highly explosive astrophysical phenomenon 
known as the Quark Nova (QN) pE]. Calculations have revealed that the energy 
released could be as high as 10 53 ergs |0]. Such excesses of energy could account 
for the mechanism behind poorly understood phenomena such as super-luminous 
supernovae J5J. Furthermore, the energetic ejecta of such events could trigger nuclear 
reactions leading to heavy element nucleo-synthesis [6], p-process element synthesis 
[Zj, and spallation processes 0[9]. Thus proving that the transition from a NS to a 
quark star is an explosive one will have transformative impact to astrophysics. 

However, making use of QCD concepts in the context of compact objects faces 
challenges. For one, the thermodynamic regimes appropriate for QCD in astrophysics 
are non-perturbative. For example, equations of states (EoSs) of nuclear/quark mat¬ 
ter that include QCD physics in studies of compact stars cannot be constructed from 
first principles. Instead phenomenological EoSs are used. Although these equations 
have been relatively successful in modelling certain aspects of compact stars, different 
EoSs models have important distinctions that would manifest in different dynamics 
in astrophysics, in particular in the interior of NSs. 

How do we choose amongst all the possible phenomenological equations of state? 
Perhaps astrophysics might be of help. Astrophysical observations that correlate 
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with specific EoS models could constraint the correct equations of state. This in 
turn, would have consequences not only in astrophysics but in particle physics in 
general, for it would shed light onto the true nature of the non-perturbative regime 
of QCD. 

Perhaps, the first step in finding out the true equation of state would be in simu¬ 
lating the micro-physics behind the phase transition of a NS into a quark star. Such 
a simulation would cover the parameter space related to disparate EoS and differ¬ 
ent initial conditions. In doing so, one could accomplish two things: (i) simulate 
astrophysical signatures that could help discriminate between the parameters, and 
ultimately the equations of state, and (ii) shed light on whether a QCD phase tran¬ 
sition inside a NS could reach an explosive regime or the phase transition simmers 
quietly and slowly. 

Here we describe the Burn-UD code which has been built and used by the Quark- 
Nova group to investigate the physics of this transition. Burn-UD nm is a hydro- 
dynamic code that simulates the micro-physics of the hadronic-quark matter phase 
transition. The work in general explores the effects of different EoSs states on the 
transition, the importance of the initial conditions (prior to the transition), and the 
hydrodynamics behind the combustion process. Other points are organized as fol¬ 
lows: Section [2] lays out the assumptions used in our model. Section 3 is a very brief 
introduction to the Burn-UD code/software which is presented as a unifying platform 
for the QCD community exploring the phases of Quark Matter and astrophysicists 
working on Compact Stars. Some preliminary results of simulations run using Burn- 
UD are discussed and examined in Section [4} Section 5 looks at future work and 
finally we presents our conclusions in Section 6. 


2 Model Assumptions 

Before we describe Burn-UD, we describe our assumptions concerning the the stability 
of uds matter, introduction of the s-quark seeding, the fluidity of quarks, the reaction 
rates of the weak interaction, and the EoSs used. 

2.1 Strange Matter Hypothesis 

The strange matter hypothesis states that matter made of up, down and strange 
quarks could be the most stable state of matter [3j. If true, then it would be interesting 
to investigate the feasibility of the phase transition from hadronic to uds matter in 
the universe. If this transition is energetic and if it leads to the formation of stable 
uds stars, the implications to astrophysics are undeniable. Here, we assume that the 
strange matter hypothesis is true. This is dynamically important because it implies 
that once a critical amount of strange quarks are seeded in the core of a NS, the whole 
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system will try to lower its energy by converting to the more stable strange matter. 


2.2 Seeding of Strange Matter 

There are numerous theories of how a s-quark seed could be introduced to the core of 
a NS. Seeding has to appear in such a way that a critical number of strange quarks 
are created instantaneously. In our model, we assume an initial seed. The cause of 
the seeding itself is left as an open question. Numerous studies have pushed forward 
different seeding models, such as clustering of hyperons in the core, and neutrino 
sparking m Another model explored for seeding is the creation of strange quarks 
through DM self-annihilation [15j [16]. This seeding mechanism is further discussed 
in Section [5711 

2.3 Asymptotic Freedom and the Fluidity of Quarks 

Because our model instantiates the combustion process in the core of the NS, we 
already assume a critical density for asymptotic freedom, i.e. hadrons in the core 
have deconhned into a fluid made of quarks. We therefore model the combustion of a 
fluid consisting (ud) matter into a fluid made of (uds) matter. Moreover, the timescale 
of the strong force, which governs quark confinement, is much shorter than the weak 
interaction, which is responsible of turning ud matter into uds matter. Thus, it is 
reasonable to assume that hadrons at the vicinity of the flame have already dissolved 
into their constituent up and down quarks [lO] . 

The conversion process of ud into uds matter is modelled as a combustion process; 
the ud matter acts as a fuel that is burnt into a uds ash. However, usually combustion 
studies usually deal with a pre-mixed fuel, i.e. the fuel is mixed with the oxidant. 
In such a pre-mixed combustion, the speed of the flame is dictated by the thermal 
conductivity, and the hydrodynamics (diffusion and advection) become unimportant 
compared to the conduction of heat across the medium. 

However, in the case of the combustion of ud matter into uds matter, the oxidant 
(strange quarks instantiated in the core of the NS) is not pre-mixed. Thus, combustion 
is only activated when the strange matter oxidant diffuses into the ud fuel. When 
the strange quarks diffuse into the ud fuel, the up, down and strange quark mix will 
attempt to achieve chemical equilibrium, converting up and down quarks into strange 
quarks, releasing heat in the process. 

Because the dynamics of combustion are dominated by a hydrodynamic process, 
i.e. diffusion of the strange quark oxidant into ud fuel, and not by heat transfer, it is 
necessary to solve numerically the full system of hydrodynamic equations. We solve 
the 1-D hydrodynamical equations: 

K = -VF(U) + S(U). (1) 
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where WF(U ) are the advection components, and S(U ) the source components. 
We refer the interested reader to [TO] for the details of the hydrodynamics of the ud 
to uds combustion. 

2.4 Reactions 

Once a strange matter seed is introduced, the following reactions will try to achieve 
equilibrium. This creates a net surplus of strange quarks, resulting in the conversion 
of the ud fluid into uds matter: 


d —y u + e~ + v e 

(2) 

u + e~ —> d + u e 

(3) 

s —y u + e~ + v e 

(4) 

u + eP —y s + W e 

(5) 

u T d —>• u + s 

(6) 


We take the rates for these reactions from See also m- 

2.5 Neutrinos 

The uds mixture will attempt to reach equilibrium through the weak interaction 
reactions. As such, the reactions will release neutrinos which act as an energy sink 
that lowers the temperature behind the burning interface (see reactions above). The 
sink is parameterized and is proportional to the neutrino mean free path. Se also [TO] . 

2.6 Equations of State and Color Superconductivity 

We assume that both the ud and uds matter near the burning interface are deconhned. 
In its early version, Burn-UD utilizes the MIT bag model m for both sides of the 
interface. The EoS used is: 



l l 


h = 4 (P + B) 


( 7 ) 

( 8 ) 


Where P is pressure, T is temperature, /i* is the chemical potential for quark 
specie i, h is enthalpy, and B is the Bag constant. 

In addition, as the uds phase (the ashes in our simulation) cools we assume that 
it becomes superconducting and enters a Color-Flavor-Locked (CFL) phase [T9]. We 
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model the CFL phase through the MIT Bag Model EoS with the following prescription 

m ' 



(9) 

A 0 = 100 MeV 

(10) 

T c = 2 1/3 x 0.57A 0 

(11) 

A = 2~ 1 / 3 A 0 i/l — (T/T c ) 2 

(12) 

P-P + 3A V 

± ~ -Tnon-CFL ~r 9 

7 T z 

(13) 

Q A 2,,2 

h = 4(P + B) 2x f 

(14) 


where T, P, h, are temperature, pressure, enthalpy. T c is the critical temperature, 
where T < T c is the temperature for phase transition to CFL. A is the CFL energy 
gap, and A 0 is the energy gap for CFL at T — 0. P n0 n-CFL is defined in equation ([7]) . 

Thus, the CFL state stiffens the EoS by adding a gap term f , and therefore 
effectively lowers the energy density of the system. This in turn releases energy that 
energizes the burning front. 

Previous studies, e.g. ra. usually model CFL in equilibrium requiring /z s = /z u = 
fid for the massless quark case. However we model the moment of the phase transition 
as a non-equilibrated system where fi s < n u . Thus we add the following condition 
that dictates when the non-CFL to CFL transition happens: 


hs ^ ^ ) 

T <T r . 


(15) 

(16) 


Thus, in a non-equilibrium situation, 3A ^ cannot simply be absorbed into an 
effective bag constant, because the CFL transition only appears when the conditions 


in equations (15) and (16) are met. 

CFL is triggered as follows: first the instantiated strange quark seed interacts 
with the ud fuel, raising fi s but lowering fi u and /x<j. After some time at a point near 


the interface, the chemical potentials will almost equal each other, and condition (15) 
is met. At this point the EoS stiffens because the component 3A 2 ;/ is added to the 
pressure equation, which creates a large pressure differential across the interface that 
in turn, triggers an advection wave travelling from the uds ash to the ud fuel. 


3 Burn-UD 

Burn-UD is a hydrodynamic combustion code used to model the phase transition of 
hadronic to quark matter. The most recent version of the code is written in java (the 
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original was written in Fortran 95, [10]) and will be publicly available soon on all 
major platforms (windows, linux, Macintosh). 

Beyond the micro-physics modelled by Burn-UD that we described in the preced¬ 
ing sections, what makes Burn-UD distinct is that it has a user-friendly graphical user 
interface (GUI). The user can modify in a transparent and user-friendly way many of 
the initial conditions, including initial temperature, distribution of the strange seed 
and initial total density. Fig. 1 shows a screenshot of the Burn-UD GUI in its most 
recent version. Furthermore, Burn-UD’s GUI allows the user to change the EoS and 
turn on and off specific hydrodynamical effects such as advection or diffusion. 

Burn-UD also contains graphic modules that plot in real time the different vari¬ 
ables of the hydrodynamic simulation. The user can select to store the data for 
specific time intervals in order to seamlessly create an animation of the results. It 
also gives the user the ability to store the data in ascii hies. 


4 Results 

4.1 Preliminary Results 

Early results using the first version of Burn-UD [bO] showed: (i) typical speeds for 
the burning front are 0.002c - 0.04c where c is the speed of light, (ii) Cooling due 
to escaped neutrinos slows down the burning front, and in some cases can halt it 
completely. 

These results hint towards an explosive regime. The details are explored in [10] but 
can be summarized as follow: Neutrino cooling slows down the burning front unevenly 
across a multidimensional surface, which wrinkles the interface. This wrinkling will 
be unstable, forming a deleptonization instability. This instability could either: (i) 
collapse the hadronic layers into the quark core and thus trigger a core collapse 
explosion or (ii) accelerate the burning front further and cause it to achieve supersonic 
speeds which would lead the system into a defragation to detonation regime, triggering 
an explosion. 

4.2 The CFL and Burning Speeds 

For the case of quark matter at a particular temperature and chemical potential, the 
CFL EoS has a lower energy density than the non-superconducting Bag equation of 
state used in our previous studies. Thus by simple energy conservation considerations, 
the transition of non-superconducting uds matter to CFL matter would release energy. 
This energy excess could also aid in leading the phase transition to an explosion. 
Another important aspect of a CFL phase transition that could accelerate the burning 
is the release of a photon fireball [IT] . Photons have momentum and thus, they will 
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Figure 1: Screenshots of the Burn-UD GUI. Upper window is the dashboard and 
the panel in the left is where the input parameters can be changed. Lower window 
is the graphing module where the chemical potentials of the strange, up, and down 
quarks are being plotted. 
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Figure 2: Position vs Time of the flame for a non-CFL (eq. (Jt]) ) EoS and a CFL 
(eq. (13)) EoS. Position is measured as the distance of the flame from the center of 
the NS. The CFL flame is much faster than the non-CFL case. 


transfer momentum to the external hadronic layers that would in turn accelerate the 
burning. 

Several tests for the CFL case have been run, and the results seem promising. 
Figure [2] shows the position of the burning front versus time in nano-seconds. CFL 
stiffens the EoS which in turn creates an advection wave moving towards the unburnt 
ud fuel, which accelerates the burning front at even higher speeds. The increase in 
speed at time ~ 3 nano-seconds in Figure [2] occurs when the condition (fi u — /j s ) < A 
(see eq. (15)) is met. Furthermore, in a multidimensional situation, the onset of 
CFL transition could further exacerbate the deleptonization wrinkling instability by 
further increasing the pressure gradient across the interface, by the simple fact that 
a CFL transition would increase the pressure in the cooler regions of the interface. 
Tests for the CFL case are ongoing to verify these claims. 
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5 Future Considerations 


5.1 Dark Matter Annihilation and Strange Quark Seeding 

A possible source of the strange quark seed may be dark matter (DM) annihilation 
mm- The distribution of DM in galaxies allows for it to be accreted gravitation¬ 
ally onto a dense and massive astrophysical object such as a NS [IT]. Figure [3] shows 
a schematic diagram representing the stages of DM accretion and self-annihilation 
to produce strange quarks. One scenario explores the self-annihilation of WIMPs 
(weakly interacting massive particles) inside NSs. They inject the NS with energetic 
photons that in turn could deconfine hadrons into bubbles of ud quark plasma. Then 
non-leptonic weak processes could create strange quarks from these bubbles, giving 
rise to uds seeds m flfi] . Whether the amount of s-quark seeding is enough (i.e. 
reaches a critical mass) to trigger the hadronic-to-quark-matter phase transition re¬ 
mains to be confirmed making Burn-UD an ideal tool/software for such investigations. 

Burn-UD is flexible enough to be used as a tool to study the strange quark seeding 
produced by DM annihilation. Burn-UD allows the user to input their own function 
for strange matter initial distribution. Different seeding models might imply different 
functions for the initial distribution of the strange matter seed. The fact that Burn- 
UD allows the user to input their own strange matter seed initial distribution makes 
Burn-UD useful for the study of many seeding models, including models for DM. 

5.2 Mass of the Strange Quark 

Up and down quarks have similar masses, yet the strange quark has a mass that is an 
order of magnitude higher. This difference between the up, down and strange mass 
should have consequences in the dynamics of the combustion. At the moment, the 
EoS used for Burn-UD assumes the quarks are massless. The only dependency in 
Burn-UD to the quark mass lies in the rates mentioned in section [4], where a strange 
quark mass of 150 MeV only affects the results very weakly as explored in [lOj. 

Different masses for quarks would affect the binding energy necessary for the 
assembling of quarks. Therefore the creation of massive quarks would deplete the 
burning front of energy and slow down the burning speed, ft is important to include 
the quark masses in future studies to investigate how the results of the Burn-UD 
simulations are affected. 

5.3 More Sophisticated Equations of State 

Burn-UD models the flame micro-physics by imposing the MIT Bag EoS on both 
sides of the interface, i.e. for both the uds ash and the ud fuel. Thus the EoS at 
both sides of the interface shares the B parameter as well. This means that at the 
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I: Energy loss 



II: Capture 

Q Dark matter 
O S-quarks 
O Nucleons 



Figure 3: Four main stages of DM capture by NSs. (I) Shows the DM energy loss by 
scattering with nucleons at high density inside the core of the NSs. (II) After losing 
energy, the DM particles are accreted in the core of the star and can’t escape the 
gravitational pull of the star. (Ill) The density of DM has increased enough so that 
self-annihilation of DM happens in the core of the star. One of the possible products 
of this self-annihilation is the strange quark. Finally in (IV) the density of strange 
quarks increases, effectively seeding the hadronic-to-quark-matter phase transition 
process. 
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region near the interface, where quarks have already deconkned, the behaviour of the 
fluid becomes independent of the B parameter. The speed of the burning front then 
becomes a function of the binding energy released during the creation and destruction 
of quarks. 

The MIT Bag Model is very popular because despite its simplicity, it manages to 
successfully reproduce some of the properties of hadrons. It also models asymptotic 
freedom in an intuitive way making it useful for simulating the properties of bulk 
quark matter. Yet such simplifications might miss some important dynamics. One 
of them is chiral symmetry breaking. A popular model that includes chirality is 
the Nambu-Jonas-Lasinio (NJL) model for QCD. One important difference in the 
NJL picture, is that strange quark matter is not the ground state of all matter [18] . 
However, the NJL model, like the MIT Bag Model, is phenomenological and as such 
its predictions are not completely conclusive either. In future releases of Burn-UD, 
the code will incorporate several EoS’s, including NJL, to compare their effects on 
the simulation. 

Finally, a more realistic simulation would include a mixed phase EoS as well. At 
the moment, deconkned quarks that act as a kuid are only considered. However, 
the dissolution of hadrons into quarks, if it is a krst order phase transition, would 
change the pressure and therefore preshock the ud fuel to a certain degree, which 
would change the dynamics of the combustion simulation. 

5.4 More Sophisticated Neutrino Transport 

As shown in |2.5[ Burn-UD currently models neutrino transport as a parameterized 
energy sink that saps the matter near the interface of energy. Although this model 
addresses many of the larger and more qualitative aspects of neutrino transport, a 
more accurate model would actually implement the time dependent Boltzmann equa¬ 
tions of neutrino transport. In the Boltzmann equations, neutrino kux is modelled 
as a density function that evolves across time and space. 

Having an accurate model of the neutrino transport is vital because the dynamics 
of these neutrinos are directly related to the deleptonization instability occurring at 
the burning interface. A sophisticated neutrino transport model will be one of the 
deciding ingredients in the discovery of whether the deleptonization instability leads 
to supersonic dekagration to detonation burning, or core-collapse Quark-Nova. 


6 Conclusion 

We hope that Burn-UD will evolve into a platform/software to be used and shared 
by the QCD community exploring the phases of Quark Matter and astrophysicists 
working on Compact Stars. Burn-UD okers the possibility of studying the time- 
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dependent phase transition of hadronic matter to quark matter. This software models 
hydrodynamically the combustion front (i.e. the flame’s dynamics) which gives a 
physical window to diagnose whether the combustion process will simmer quietly and 
slowly, or transitions from deflagration to detonation. However, it is important to note 
that there exists solutions where the burning front is halted due to the deleptonization 
instability which make the system undergo a core-collapse QN explosion. Whether 
via a Deflagration-Detonation-Transition or core-collapse, such an energetic phase 
transition would have many consequences in high-energy astrophysics and could aid 
in our understanding of many still enigmatic astrophysical transients. Furthermore, 
having a precise understanding of the phase transition dynamics for different EoSs 
could aid further in constraining the nature of the non-perturbative regimes of QCD 
in general. 
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